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ABSTRACT 


As the target radar eross seetion (RCS) eontinuously deereases, the need for high- 
resolution high-gain radar inereases. One approaeh to high resolution is to use distributed 
subarray antennas (DSAs), because of limited surface available on many radar platforms. 

The concept of distributed subarray antennas is examined for both Multifunction 
Array Radar (MFAR) and Very High Frequency (VHF) applications. By combining 
distributed subarrays located on the available areas of a constrained platform, the MFAR 
and VHF DSA can achieve increased resolution and potential reductions in cost and 
complexity compared to a conventional array. The two-way pattern design of DSA 
effectively suppresses the undesired grating lobes by using separate transmit and receive 
antennas. By the pattern multiplication principle, the grating lobes in the receive pattern 
have been suppressed by proper null placement of subarray in the receive and transmit 
patterns. 
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I. INTRODUCTION 


A, REQUIREMENTS AND OBJECTIVES 

The platform design of future surfaee combatants has changed dramatically with 
the advent of stealth technology to reduce platform signatures. The change in design 
philosophy is evident in the proposed new DD-21, as illustrated in Figure 1. The 
traditional small-integrated superstructure of vessels allowed few areas for sensors and 
weapon systems. On the other hand, the threat from the air or surface is ever increasing 
such that the performance of shipboard radar needs to provide long-range detection and 
accurate tracking. This implies high gain and physically large radar antennas. 

Since it is difficult to find a sufficient area for a large array onboard a ship, it 
might be possible to use several smaller noncontiguous (separated) areas (subarrays) and 
then process the received signal coherently. The subarrays may be far apart in terms of 
wavelength and therefore grating lobes occur. Grating lobes are undesirable because they 
reduce the antenna efficiency, cause ambiguities in angle measurements and make the 
radar more susceptible to jammers. They also complicate clutter processing. 

The main objective of this thesis is to investigate a conceptual design of 
integrated antenna apertures, which are composed of several distributed noncontiguous 
subarrays. Together the entire set forms a Multi-function Array Radar (MFAR) antenna 
that provides multi-tasking with high-resolution high-speed data collection 
simultaneously. This thesis will show the trade-off in performance and cost for several 
possible antenna concepts. 

One possible approach to achieving a large distributed array that is free of grating 
lobes is to design separate receive and transmit antennas. For example, the transmit 
antenna might consist of distributed subarrays. Then a separate receive antenna would be 
used whose pattern has nulls in the direction of the transmit array’s grating lobes. 
Therefore the round trip or “two way” gain pattern would be free of grating lobes. 
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Figure 1. The Artist’s concept of the DD-21 [Ref. 1] 


This approach is difficult enough for narrow band radar. However, new antenna 
designs must provide wideband performance. Furthermore, due to the limited space, it is 
desired that the antenna serve as many systems as possible: several radars, 
communications systems, and electronic warfare systems. Ideally, the antenna should 
have the following characteristics: 

1. Capability to execute different tasks in rapid sequence 

All sensors and electronic devices need to be integrated into the small topside of 
the platform, and the integrated MFAR should have surveillance, tracking, identification, 
lire control, missile guidance and communication capabilities. The bandwidth of the 
MFAR, therefore, must be wide enough to provide good range resolution and to satisfy 
the requirements of different systems. The aperture size of the MFAR must be electrically 
large enough to give fine angle resolution and high directivity. The side lobes of the 
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antenna pattern, of course, have to be low enough to help distinguish the target from the 
natural or man-made clutter. 

2. Capability for adaptive digital beamforming (ADBF) 

Digital beamforming (which requires a modulator or demodulator at each 
element) is preferred. However, the large number of array elements required for high 
resolution reduces the possibilities of element-level digital beamforming considerably. 
The subarray-level ADBF could provide multiple beams and electronic scanning to carry 
out the multi-tasks simultaneously using a minimum number of control elements. The 
MFAR tracks targets accurately with low sidelobe sum and difference beams, which are 
required for monopulse tracking. There are also some other important features that could 
be provided by ADBF: 

(1) Self-calibration and error correction 

(2) Adaptive nulling of unwanted interference or jamming signals 

(3) Spatial control of the radiated energy for Low Probability of Intercept 
(LPI) operation 

3. Affordability with acceptable performance 

The purpose of Multifunction Subarray Radar (MFSAR) is to fulfill a variety of 
different tasks. Generally the cost is proportional to the number of radiating and control 
elements, as well the complexity of construction and maintenance. For a modern high 
performance radar antenna, the basic acceptable antenna parameters are the following 
[Ref 2]: 

Antenna gains: 30 dB (minimum) 

Target tracking accuracy: 5 mrad 
Azimuth beam steering: ±45° (Azimuth) 

Azimuth/elevation coverage: 3607 60° 

Peak sidelobe level: -45 dB 
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B, PREVIOUS APPROACHES 

There are two basie approaehes that have been eonsidered as a solution for a 
wideband multifunction system from the standpoint of the array aperture. The first is a 
segmented approach, where each segment or component is dedicated to a separate system 
and bandwidth. A segmented aperture provides better performance for each individual 
function by separating the size and position of the antenna. The second approach, a 
shared aperture, however, provides multiple frequency operation in the same aperture. 
Generally the performance at a single frequency is not optimum. 

The proliferation of advanced sensor and communication systems onboard 
military platforms has led to a multitude of systems. A U.S Navy Aegis cruiser, for 
instance, has over 100 antennas and the number is expected to rise as new systems are 
added [Ref 3]. A reduction of the number of antennas is possible using the shared 
aperture concept. 

Ideally a shared aperture of the type used on ships should operate over a wide 
bandwidth, for example, 10 MHz to 10 GHz (three decades). If the aperture is shared by 
separate narrow band systems, then the antenna does not have to operate over the 
complete continuous range of frequencies, but only at the operating “sub-bands” of the 
component systems. These types of antennas are referred to as multi-band. An example is 
the Multifunction Electronically Scanned Adaptive Radar (MESAR) [Ref. 4] that 
employs a dual frequency antenna. The work in MESAR MEAR programs, which began 
in 1977, showed significant advantages for a radar operating at dual optimal frequencies. 
This study led to the selection of 1 and 10 GHZ for operating frequencies used for 
surveillance and tracking functions. The aperture is comprised of two sets of antenna 
elements, an open-ended waveguide and a dipole, for the low and high frequency bands, 
respectively. The proposed system will therefore perform these important radar functions 
at maximum efficiency. 
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C. MULTIFUNCTION SUBARRAY RADAR 

Modern phased array radars are used in a multifunction fashion, with the required 
functions being search, tracking and fire control. A major impetus of using MFSAR is 
cost reduction. For example, if phase control is at the subarray level and each subarray is 
formed by 5x5 elements, then the total number of transmit/receive modules is reduced by 
a factor of 25 for contiguous subarrays. With widely spaced subarrays, not is there a 
savings in control elements, but also higher resolution. When the subarrays are widely 
distributed the potential resolution increases significantly. The disadvantage is the 
occurrence of grating lobes if the subarrays are widely spaced compared to the 
wavelength. The features of MFSAR are summarized in the following sections. 

1. Advantages of MFSAR 

The major advantages of using multifunction subarray radar are the following; 

• More efficient scheduling of the track and search functions compared to 
separate systems (track-while-scan, TWS capability). 

• Rapid steering of the beam to the desired direction as needed. 

• Formation of search and track beams with maximum flexibility. 

• Savings in the number of control elements and decrease in the complexity 
of digital beamforming. 

• Increase in the angular resolution by distributing subarrays along the 
superstructure of vessel. 

2, Disadvantages of MFSAR 

The major disadvantages of using multifunction subarray radar are the following: 

• Compromise in the performance of individual functions to obtain an 
optimum balance in the component tasks. 

• Extra grating lobes introduced by large subarray spacing. 

• Low beam efficiency from the presence of grating lobes. 
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D, SCOPE AND ORGANIZATION 

1. Scope 

In this research, some methods of correcting the above disadvantages are 
examined. The grating lobe problem is severe for subarray spacing over a wavelength. 
However, it is possible to suppress the grating lobes through control of the subarray 
factor. In other words, the nulls of the subarray pattern are placed at angles where the 
grating lobes occur. There is a limit to the effectiveness of this method, as the gaps 
between subarrays get larger. Some additional improvement is achieved if separate 
transmit and receive antennas are used, and grating lobes are allowed for only one of the 
two antennas. For example, if the transmit antenna has grating lobes, then the receive 
antenna will not. Furthermore, the nulls of the receive antenna can be placed at transmit 
antenna grating lobe locations. Consequently the two-way gain pattern will have no 
grating lobes. 

An additional improvement may be possible with digital arrays that have 
amplitude control at each element. The transmit array factor and subarray factor, as well 
as the receive array, can be phase and amplitude weighted to give low sidelobe 
performance. 

This thesis demonstrates that distributed antennas of MFSAR can provide better 
performance with fewer elements than the conventional MFAR by using separated 
subarrays. 

2, Primary Research Question 

There are two related research questions addressed in this research; 

a) how to increase the angular resolution of shipboard radar through wideband 
distributed subarray antennas, and 

b) how to suppress the grating lobes in the pattern of widely spaced subarrays by design 
of the subarray pattern or design of a two-way pattern. 
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3. Organization 

Chapter II provides an overview of the array theory that is used in this thesis. A 
definition of the array factor (AF) is given. Antenna parameters such as beamwidth, 
directivity, aperture efficiency, grating lobes and pattern scanning, are presented. The 
assumptions and limitations are introduced. 

After a description of various subarray configurations. Chapter III discusses the 
different methods used in the conceptual design of an MFSAR antenna and the concept of 
digital beamforming will be discussed in more detail. Some ship design examples are also 
presented. 

Chapter IV provides a summary and conclusions, followed by suggestions for 
further research into distributed subarray design. 

Appendix A provides a glossary of terms and abbreviations used throughout this 
thesis, and Appendix B provides a listing of MATLAB codes used in the pattern 
calculations. 
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II. SUMMARY OF ARRAY THEORY 


This chapter discusses the fundamental theory of phased arrays. The material can 
be found in most books on antennas, such as [Ref. 5] and [Ref. 6]. 

A, ARRAY FACTOR (AF) 

An array is a collection of smaller, usually identical antenna elements that are 
excited with complex voltages or currents to obtain a desired radiation pattern. For most 
applications the elements are arranged in a periodic grid; for example, in two dimensions 
it could be a rectangular lattice. The array can be divided into smaller groups of elements 
called subarrays as shown in Figure 2. For simplicity, only a linear array is shown, and 
the array elements are assumed isotropic. 


i!i=90‘ 


< 


I 


Isotropic element 



^= 0 * 


Figure 2. Example of an array of contiguous subarrays, M = 5 , N = 3. 
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The radiation pattern of an antenna array is the veetor sum of the electrie field 
intensity all antenna elements. The array faetor for uniform subarrays, eaeh eomposed 
of M isotropie elements along the x -axis is given by 


V-l M-1 


AF(e.^) = yyi,J^e 


2)r . In 

jn^^D^{u-u^) jm — d^{u-u^) 


«=0 m =0 


M-\ . 2;r , . . 

7 J e 

/ j m 


S 

II 

o 

n=0 


( 1 ) 


= AFXd,(l>)xAF^{d,(l>), 

where 

N,M = number of subarrays, number of elements in a subarray 
l„,J subarray, element weights (excitation currents or voltages) 
D^,d^ = subarray center-to-center spacing, element spacing 
u= sin 0 cos ^, direction cosine in spherical coordinates 
= sin 0^ cos {0 ^,is the scan direction) 


AF^ = subarray pattern determined by the arrangement of elements in a subarray 
AF^ = configuration pattern determined by the arrangement of subarrays. 

Equation (1) assumes that all subarrays are identical and the elements have 
uniform spacing. The pattern of the array is then the product of the element factor and the 
array factor. Isotropic elements have been assumed and, therefore, the element factor 
(EF) is 1. The use of subarrays allows the array factor to be separated based on the 
geometry and excitation of subarrays. 
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B, ANTENNA PARAMETERS 

1. Beamwidth Between First Nulls (BWFN) 

The beamwidth between the first nulls is defined by the subtended angle of the 
mainlobe in the antenna pattern. If the element and subarray weights in Equation (1) are 
assumed unity, then the normalized antenna pattern is given by: 


where 






M sin(^) N sin( y) 




Itt 

T 


d^{u-uj 


( 2 ) 


y/ = —D^{u-u^). 

The subseript “norm” denotes the normalized array faetor. The BWFN is determined by 
the first zero of the numerator. Sinee the phase term of the subarray is changing much 
faster than the phase of elements ( ND^ » Md^ ), the BWFN is determined mainly by the 

phase of subarrays given by u-u^ =±1/A®^ so = 2sin“'(l/in the plane 
= 0 . The BWFN, therefore, is determined by the entire aperture size (i.e., the distance 
between the edges of the two farthest separated subarrays, L ~ ND^ ). 

2, Half Power Beamwidth (HPBW) 

The half-power beamwidth of a linear array is defined as the angular separation 
between two points, one on each side of the main beam maximum, at which the power is 
reduced by half The HPBW is determined mainly by the subarray configuration pattern. 
For a half-wavelength spaced linear array with uniform excitation, the HPBW can be 
approximated by 0hpb„, ~ 50.8/1/ L in degrees, where L ~ ND^ is the entire length of the 
aperture. 
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3, Directivity and Gain 

Directivity is defined as the ratio of the maximum radiation intensity in the main 
beam to the average radiation intensity [Ref. 7]. For a two-dimensional uniform 
rectangular array of isotropic elements, the directivity can be approximated by 
D = A7r7]L^Ly cos 6^ / 2} , where 7] is the aperture efficiency, L^,Ly are the dimensions of 

the array in thex,y directions, and 0^ is the scan direction. This approximation is good as 
long as there are no grating lobes in the visible region. Subsequently it is assumed that 
there are no other losses other than those due to amplitude tapering. In this case the 
directivity and gain are equal. 

4. Aperture Efficiency 

Aperture efficiency is a measure of how efficiently the antenna physical area is 
utilized. If the element or subarray amplitudes are not equal, the aperture efficiency is 
given by the taper efficiency. For a planar array in the x- y plane 


71 = 






^.Jv,ZZK. 


n=\ m=l 


( 3 ) 


where and are the number of elements in the x and y directions, respectively, and 
is the amplitude of the mn th element weight. The weighting can be applied to either 
the AF^ or AF^ summation. 

5, Grating Lobes 

Grating lobes occur when more than one period of the array factor appears in the 
visible region (±90°); when either the element or subarray spacing in the Equation (1) is 
more than one wavelength. For MFAR it is likely the subarray spacing is larger than one 
wavelength, and thus grating lobes will exist, as shown in Figure 3. In theory, the grating 
lobes could be suppressed by using either unequally sized or randomly spaced subarrays, 
at the expense of higher complexity, average sidelobe level increase, and difficulties in 
applying ADBF [Ref 8]. Another possible means of suppressing grating lobes is to use 
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two separate antennas for transmit and receive (i.e., a quasi-monostatic radar). The two- 
way gain occurs in the radar range equation for the signal-to-noise ratio {SNR) 


SNR = 


Pr_ P,Gfi,Xc7 
N {ATtfR^KT^B^ ’ 


( 4 ) 


where is the received target signal, is the transmit power, a the target RCS, A, the 
wavelength and R the range. The product KT^B^ is the noise power, K is Boltzmann 
constant, system noise temperature and B^ radar bandwidth. Grating lobes can be 

allowed in the receive pattern and then eliminated by placing transmit nulls in those 
directions. The cost is usually a reduction in transmit antenna efficiency. 



Figure 3. Grating lobes produced by subarrays = 5A = 0.5/1, = 5,M = 5). 
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6, Electronic Scanning 

Electronic scanning is defined as a method of positioning a beam with the antenna 
aperture remaining fixed. The basic electronic scanning techniques are phase shifting, 
true time delay, frequency scanning, and feed switching. Some possible control 
configurations for beamforming are shown in Figures 4 through 6. To avoid the problems 
of beam squinting and broadening that occur over wide frequency bands, it is desirable to 
have true time delay at each element. This is difficult to achieve using conventional 
microwave beamforming. However, it has been accomplished using photonic devices 
(optical fiber) and through signal processing in digital antennas. Subsequently, true time 
delay is assumed at each element. 



PHASE 

SHIFTERS 


SUBARRAYS 

IN/OUT 


TRUE TIME 
DELAY 


Figure 4. True time delay at the subarray level with phase shifters at the element 
level. 
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Figure 6. Digital beamforming. 
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c. 


RADAR SYSTEM CONSIDERATIONS 


The radar range equation (4) gives the fundamental relationship between radar 
design parameters. A main objective of this thesis is to propose a possible conceptual 
design of distributed subarray antennas, using separate transmitter and receiver antennas, 
which is a form of quasi-monostatic radar. 

1. Maximum Detection Range 

From Equation (4), the maximum detection range of the radar for a given signal- 
to-noise ratio can be rewritten as 

^ 4 _ _ P,G,A^C7 

- {A7rf{SNR)^KT^B^ {A7rf KT^BSSNR\ 

{47rfKT^.FSBT)f^{SNR)„’ 

where 

Pgvg - average transmitter power 
= receiver effective aperture area 


= receiver noise figure 

B = receiver bandwidth 

T= pulse width 

fp = pulse repetition frequency 

(SNR)^ = minimum signal-to-noise ratio after n -pulses are integrated. 

This equation illustrates several important tradeoffs in the design of radar 

(1) Power-aperture product: One of the most important measures of the 
capability of long-range surveillance radar is the product of average power 
and effective aperture size (F^^gA^J , which controls the amount of power 

transmitted by the aperture. If the aperture size could be larger, the average 
power required for detecting the target at a given range would be decreased 
when all other factors are constant. 
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(2) Frequency dependence; Although the frequency does not appear explicitly, it 
is easier to achieve high power at low frequencies because high voltages are 
applied to larger distances and breakdown can be avoided. Furthermore, 
atmospheric loss is less at lower frequencies. 

2. Angular and Range Accuracy 

Range and angular resolution provide improved accuracy, since accuracy is 
inversely proportional to signal bandwidth and directly proportional to beamwidth as 
[Ref 9] 




e, 


HPBW 


KJ2SNR, 


( 6 ) 


where 


<7^= angular accuracy 


Ohpbw ^ antenna beamwidth 


= gradient of the difference beam in the monopulse antenna configuration. 

Once again, the angular accuracy is proportional to the beamwidth and inversely 

proportional to square root of SNR in this equation. 

Since the range resolution is AR = ct/2 = c/2B [Ref. 3], the range accuracy is 
improved by having a wide instantaneous bandwidth. 
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D, SUMMARY 

This chapter has presented the basic theory of arrays and discussed some 
fundamental radar system design tradeoffs. This researeh is primarily coneemed with 
increasing radar resolution by an increase in the antenna aperture. Because large 
unobstructed smooth surfaces are limited on ships and aircraft, it is necessary to construct 
large arrays by combining signals from distributed subarrays. Unfortunately, this results 
in grating lobes. 

In order to suppress the grating lobes the subarray sizes, loeations and excitations 
can be adjusted. Furthermore, separate transmit and reeeive antennas can be employed, 
and designed so that the two-way pattern has suppressed grating lobes. These last two 
approaehes are examined in subsequent chapters. To simplify the analysis the following 
assumptions are made: 

1. The element faetor is negleeted, 

2. Generally, a linear array is used with x being the array axis, 

3. The antenna is at the center of a spherical coordinate system, where the 
x-y plane is the earth’s surface and z the zenith direction, 

4. Mutual coupling is neglected, and 

5. A eonstant frequeney or time-harmonie wave is assumed (e^"^Uime 
dependence), and therefore phaser quantities appear in the equations. 
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III. DISTRIBUTED SUBARRAY ANTENNA 


A. SUBARRAY METHODS AND CONFIGURATIONS 

The concept of subarraying arises from the requirement of modem radar for high 
resolution. When the targets to be tracked have very small RCS, and the radar must have 
fine angle resolution to track multiple closely spaced targets, the antenna apertures have 
to be electrically large. i Because of the large number of array elements required, 
grouping of elements to reduce the beamforming complexity and control cost is 
inevitable. The subarrays can be divided into several types as shown in Figure 7. They 
are discussed in the following section. 

1. Subarray Configurations 

In terms of control there are generally four types of subarray configurations, as 
shown in Figure 7: 

1) Amplitude and phase control at each element (Figure 7(a)) — this allows 
complete control of both AF^ and AF ^, 

2) Phase only control at each element (Figure 7(b)) — this allows scanning of 
both AF^ and AF ^, 

3) Amplitude and phase control of the subarray pattern AF^ (Figure 7(c)) and 

4) Amplitude and phase control of the array factor AF^ (Figure 7(d)). 


1 Synthetic aperture radar (SAR) is another solution to the resolution problem, but it has its own 
disadvantages and limitations [Ref. 10]. 
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AMPLITUDE/PHASE (a) 


PHASE-ONLY (b) 



SELECT ELEMENTS (c) SUfiARRAY WEIGHTING (d) 

Figure 7. Typical subarray control configuration [From Ref. 11] 

Amplitude and phase control at each element is undoubtedly the most costly but 
also the most desired in terms of beam pointing and sidelobe level (SLL). The second 
alternative is to control element phase only, which is very attractive since in a 
conventional phased array the required controls are available at no extra cost. However, 
the grating lobes in the DSA approach severely distort the pattern, and a phase-only 
correction is not very effective. In this thesis, it is assumed that complex weights can be 
applied at both the element level and the subarray level. Furthermore, the weights have 
true time delay behavior.2 


2 As a practical matter, if complex weights can be applied at the element level, then there is no need 
for weighting at the subarray level. 
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2, Basic Properties of DSA Pattern 

Due to the large spaeing between subarrays in a periodic DSA, there are many 
grating lobes in the visible region. For example, if the spacing between subarrays is 5 
wavelengths, there will be 10 grating lobes in the [-;r/2,;r/2] visible region, that are 

located at the angles sm“'(±/»/5) radians (/> = 1,2,3,---) as shown in Figure 3. The 
following are some basic characteristics of the periodic subarray configuration: 

• The array pattern of identical periodic subarrays can be separated as a 
multiplication of a subarray factor (AFJ and the subarray geometry and 

excitation factor (AFJ as shown previously in Equation (1), 

• The number and the intensity level of sidelobes between two adjacent 
grating lobes are dependent on the number and position of the subarrays. 
Similar to the uniform linear array, there will be A - 2 sidelobes with a 
peak intensity level about -13 dB in the unity weighting case, and 

• The number of subarrays does not have any effect on the number of the 
grating lobes as long as subarrays are spaced approximately equal. The 
intensity level of the sidelobes is decreased with an increase in the number 
of subarrays. This implies that having multiple small subarrays that are 
spaced uniformly gives better performance in both SEE and grating lobes 
than having two large subarrays spaced far apart. 

3, Methods for DSA Pattern Design 

There are two approaches to reducing grating lobes. The first is based on pattern 
multiplication, as illustrated in Equation (1). The grating lobes in one factor can be 
suppressed by placing nulls of the other factor coincident with the grating lobes. The key 
to this approach is that the grating lobes and nulls are periodic in arrangement. A second 
approach to reducing grating lobes is to use an irregular spacing or unequal sized 
subarrays, thereby reducing the peak grating lobes by redistributing the energy into the 
sidelobe regions. 
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The process of pattern design can be divided into two stages. The first stage is to 
lower the SLL and eliminate the grating lobes as much as possible by synthesis of the 
appropriate AF^ and AF^ functions. For the subarray spacings of interest the first pair of 

grating lobes (one on either side of the main lobe) are not easy to suppress because of 
their narrow width and close proximity to the main lobe. The second step is to suppress 
the remaining grating lobes by either a high contrast (i.e., low sidelobes relative to the 
mainlobe) transmit pattern or by specific placement of the nulls in the transmit pattern. 

There are several simple methods used in the pattern synthesis of DSA in this 
thesis. The process of direct nulling by pattern synthesis is discussed next. 

a. Pattern Multiplication Principle 

Since the antenna is composed of periodic subarrays with the same 
spacing, their pattern could be predicted from Equation (1) and the null locations 
manipulated. The angular directions of grating lobes in the configuration pattern are 
given by 


u = u^± 



(V) 


where p is a non-zero integer. The nulls of the uniform subarray pattern AF^ are located 
at 


u = u^± 


pX 


( 8 ) 


The condition for suppressing grating lobes is D^/d^=M. Unfortunately this is exactly 

the condition of contiguous subarray, which contradicts our DSA approach. Since the 
nulls of the uniform subarray pattern are located in 2;r/M angular increments, if the 
element number M is chosen correctly, some ratio of grating lobes in the configuration 
pattern will be suppressed. For example, if D^ld^ = 5XI0.5X , and let M = 5 , the even 
integer numbers of grating lobes will be suppressed as shown in Figure 8 (a). 
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This only solves half of the problem. The two-way pattern synthesis gives 
additional freedom to use the nulls of the transmit pattern to suppress the remaining 
grating lobes in the reeeive pattern. The two-way pattern is defined as 

= AFt. X AF^ = AF,, X AF^ X AF ^. (9) 

For example, the reeeive pattern of 16 weighted subarrays, and transmit pattern of 20 
weighted elements, spaeed in 0.5A is shown in Figure 9. In this example the transmit 
antenna has a Chebyshev distribution. 
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Receive pattern of 16 subarrays spacing 5X 



Pattern of singie subarray which is composed by 5 eiements 



Pattern of 16 subarray configuration 



(b) 


(c) 


Figure 8. The grating lobes suppressed by 5 element subarrays pattern, (a) total 
receive array pattern, (b) subarray factor AF ^, and (c) subarray configuration factor AF^. 
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Two Way pattern of 16 subarrays 



(a) 


(b) 


(c) 


Figure 9. Example of pattern multiplication principle (a) two-way pattern , 

(b) receive pattern AFj^, (c) transmit pattern AFj.^. 
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b. Modified Minimax Algorithm to Find the Geometry of Lowest 
Sidelobes by Perturbation of Subarray Location 

There are many different approaches to the design of low sidelobe patterns 
in phased array systems. Since the goal of the first step is to minimize the sidelobes and 
grating lobes, one such approach is to displace the subarray locations based on the Dolph- 
Chebyshev theorem [Ref. 12]. 

Every subarray considered is identical (same AF^) and has symmetric 
weights about the center of the array. The AF^ can be represented as 


where 


^-1 271 D 

AFfu) = '^0Ci^ cos(- -u) , 0 < M = sin^ < 1, 

A, 


k=0 


( 10 ) 


N 

K = — when N even, 
2 


A^ + 1 
2 


when N odd, N is the number of subarrays and 


= distance between the Ath subarray and 0th subarray. 

For odd length arrays the weights to be applied on the individual subarrays are found 
from L=(Xn, f=I „=oc„l'^ and in the even case / = / = a / 2. 

By fixing the location of the center and edge subarrays, the lower bound 
on the beamwidth of the main lobe of AF^ is approximately ~ 0.886/1/ L where L is 

the length of aperture, as the initial guess. The Minimax algorithm is used by setting the 
lower and upper bound for location of other subarrays, and the obtaining subarray- 
optimized positions by minimaxing the Chebyshev basis function in the region Uq<u<\. 
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The application of this method to the 16 uniform subarrays of 5 uniform 
elements each is shown in Figure 10. Compared to the original receive pattern of 16 
Hamming-weighted subarrays with equal spacing in Figure 8(a), although the peak 
grating sidelobes are decreased from 4 dB to 11 dB, the average sidelobe level becomes 
much higher because of the perturbation in locations. Also note that the narrow 
beamwidth is similar to the uniform weighted subarray. The locations of the subarrays in 
wavelengths are Z)„= [0 2.5000 11.6208 15.1714 19.2073 24.6143 29.2136 

35.0000 40.0000 45.7864 50.3857 55.7927 59.8286 63.3792 72.5000 

75.0000]. 



Figure 10. One-way array pattern of 16 location-perturbed uniform subarrays 
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The same Minimax algorithm ean also be extended to synthesis of the 
weights of the desired transmit pattern. It will find the suitable weights of the transmit 
elements for seleetive nulling angles and the speeified sidelobe level, whieh match the 
grating angles in the receive pattern. The objective function becomes [Ref. 13] 

min{max^ |lfi(M,)(Z)(M,) -R(m,))|} = d (11) 

where 

u = ;rsin^ , the set of spatial frequencies the response is optimized 
W = the error weighting 
D = the desired response (intensity level) 
d = deviation between the pattern and desired response 

X 2 £) 

P{u,) = cos(4m) , 4 = . 

k=0 ^ 

This theorem could be formulated as a matrix equation Aa = D where 



COS * * * 

COS^j^Mq 

d 

W(u,) 

A = 

cos 

COS^^Mj 

-d 

W(u,) 


COS^qM^^j 

COS^^M^+J 

W'K.,). 

II 



D = [ 

Z)(Mo) D(u^) • 




The coefficients that form the solution to this system can then be used as the element 

weights of the transmit array. This approach was tried but the null depths were not always 
sufficient to suppress the grating lobes. 
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c. Modified Weighting Method for Distributed Subarray 

The conventional weighting distributions can be applied to nonuniformly 
spaced arrays by sampling the continuous distribution at the appropriate points that 
correspond to the element locations. The phase-shift of elements for scanning needs to be 
modified accordingly [Ref 14], 

Let the number of the subarrays be N, the number of array elements in 
the n th subarray M{n) = M (identical subarrays), the distance between the nth and 0th 
subarrays are L{n) wavelengths, and 0^ the angle from normal direction of the array. 
Then the phase-shift of the mth element in nth subarray is 

'T.tt 'T.tt 

(Pmn=L{n)—sm9,+m—d^sme^. ( 12 ) 

The modified Hamming weighting function for a nonuniform spaced subarray is written 
as 


where 


Wfjin) = 0.54 -0.46 cos[^^^^]. 


(13) 


4(n) = L{n)l nth subarray position relative to the element spacing 


=L^I d^ aperture dimension relative to the element spacing. 

A similar result could be achieved if this modified weighting method is applied to the 
elements in all subarrays, which yields a smoother sampled distribution than weighting at 
the subarray level only. However, the array factor cannot be separated because of 
different subarray factors AF^. Also 4(n) has to be changed to the element relative 

position lj{mn) in the aperture to apply the weighting function of Equation (13). 
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Figure 11 shows the reeeive pattern of modified Hamming weight 
applieation on the elements of the same subarrays. Compared to Figure 9(b), there is 
more suppression of the grating sidelobes at the angles 0 = ±23°,±53°, but almost the 
same level everywhere else. The patterns for weighted subarrays versus weighted 
elements do not differ too much in beamwidth and sidelobe level. The subarray level 
control scheme is more desired because of its simplicity. 

For reference, Figure 12 shows the pattern of 16 uniform weighted 
subarrays. The pattern of the Hamming weighted elements or subarrays is “cleaner” in 
terms of sidelobe level especially in the proximity of the mainlobe. The Hamming pattern 
has a slightly wider beamwidth. 



Figure 11. The DSA receive pattern with modified Hamming weights on array 
elements 
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Figure 12. The DSA Receive pattern with uniform weighted subarray elements 
(A = 16,M = 5,Z)^ = 5/l,J^=0.5/l) 

d. Equiripple (Parks-McClellan) Design Method of High Contrast 
Transmit Pattern 

Linear antenna arrays are in many ways analogous to one-dimensional 
digital filters. Restriction of the pattern synthesis problem to that of discrete arrays of 
finite spatial extent makes the problem similar to that of finite impulse response (FIR) 
digital filters. When the ideal time delays are used, wideband pattern synthesis reduces to 
the narrow band case, with each element’s delayed waveform receiving a single real 
weight. In classic narrow band pattern synthesis, an equiripple weighting with narrowest 
beamwidth for a given sidelobe level was proposed in [Ref 15]. 
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The Parks-McClellan algorithm is based on an iterative algorithm, which 
minimizes the maximum amplitude of the ripple (side lobes) present. By the specified 
angle of the passband (BWFN), stopband (spatial directions outside the mainlobe), 
frequency response (intensity level) and the maximum deviation from the frequency 
response the weighted least squares algorithm (WLS) provides an optimal approximation 
to the desired pattern in the least squares sense. This ensures that the power present in the 
stopband will be a minimum. The calculation of an array factor of 35 linear elements 
spaced XU using MATLAB’s built-in Remez function is shown in Figure 13, and the 
applied weights for each element are shown on Figure 14. 



Figure 13. Transmit pattern of 35 equiripple-weighted elements 
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Figure 14. Weights of the transmit elements 

B. POSSIBLE SHIPBOARD MFAR DSA DESIGN 

There are several reasons why a distributed subarray antenna might be used on a 
ship. One is the limited surfaee area available for antenna placement. At low frequencies, 
the open surface areas are small compared to the wavelength, and narrow beamwidths are 
not possible with a single contiguous array. By adding more subarrays on other areas of 
the ship, they can be processed as a DSA resulting in a half power beamwidth of 
approximately XIL radians. Any combination of transmit and receive functions could be 
used. In Figure 15, for example, 1 and 2 transmit only; 3 receives only, etc. This 
approach can be applied to frequencies where the subarray spacing is in the range of 1 to 
5 wavelengths. For an Aegis-sized cruiser this would be in the VFIF to UHF frequency 
regions. 

At higher frequencies a DSA might be used to reduce cost and weight. As 
illustrated in Figure 16 the original aperture size is large enough (in wavelengths) to 
provide a sufficiently narrow beam, XIL. Some weight can be eliminated by removing 
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selected areas of the array, as shown in Figure 17. The thinned areas between the 
subarrays can be used for other sensors; for example, radar or communication antennas at 
other frequencies. The penality is grating lobes, although they can be suppressed using 
the techniques described previously. 
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1. MFAR DSA Configuration by Thinning 

One possible antenna configuration for a shipboard DSA is shown in the Figure 
17. C-band is chosen as the operating frequency for this design. The transmit array is 
composed of the center 35x35 uniform elements (yellow), which also can be used as 
part of the receive antenna (blue squares). In theory the transmit pattern could use any 
number of elements, with the performance improving with the number of elements 
(narrower transmit beamwidth and higher gain). 

The receive antenna uses digital beamforming with 16x8 subarrays, each 
comprised of 5x5 isotropic elements. Each subarray can have any number of 
independent functions such as communication, missile guidance, sidelobe cancellation, 
etc., and they are grouped independently for each function. The following section shows 
patterns for 16 by 8 subarrays for receive and a continuous 35 by 35 element transmit 
array. 
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Figure 17. Configuration of MFAR DSA 
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2, Pattern Synthesis of Shipboard MFAR DSA 

Based on the synthesis method introdueed on Seetion A of this ehapter, Figure 18 
shows the receive pattern from coherent signal processing of 16x8 uniform subarrays. 
The grating lobes exist at 0 = 11.5°,23°,37°,53° (main beam not scanned). By 
multiplying the transmit pattern of the center 35 elements as shown on Figure 13, the 
resulting two-way pattern in the H and E-planes is given in Figures 19 and 20. The peak 
SEE is approximately ^5 dB. The grating lobes in the two-way pattern have been 
reduced to a maximum of ^5 dB below the two-way main beam gain. This demonstrates 
the effectiveness of the 2-way pattern approach. 



Eigure 18. Receive H-plane pattern of 16 Hamming weighted subarrays spaced 5A, 
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Figure 19. Two-way H-plane broadside pattern of MFAR DSA 



Figure 20. Two-way E-plane broadside pattern of MFAR DSA 
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Figures 21 and 22 show the receive patterns when the beam is scanned to 
= 10° and 60° respectively. The nulling effect of subarray factor on the configuration 
factor is unchanged with scan. The transmit pattern broadens from scanning at the same 
rate as the receive pattern broadens. Figures 23 and 24 show the two-way pattern 
scanning to 9^ =10° and 60°, Figure 25 shows the two-way pattern of multiple beams 
scanned at increments of A9 = 5°. 



Figure 21. Receive pattern of 16 subarrays when scanned to 10° from Fl-plane 
broadside 
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Figure 22. Receive pattern of 16 subarrays scanned to 60° from H-plane broadside 



Figure 23. Two-way pattern scanned to 10° from broadside 
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Figure 24. Two-way pattern scanned to 60° from broadside 



Figure 25. Multiple beams scanning with increments 5° 
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3, Shipboard HF and VHF Distributed Subarray Antennas 

The other possible applieation of DSA onboard a ship is at the lower frequeneies. 
Lower frequencies have better performance for long-range surveillance, especially on 
targets with small RCS because HF and VHF are in the resonance region for the targets 
of interest (e.g., cruise missiles). However, a major problem is how to increase the 
resolution from the limited small areas on a ship. By processing signals coherently from 
several separate subarrays, the half power beamwidth can approach XIL, where L can 
be the total length of the ship. 

Figure 26 is the side view of an AEGIS cruiser, with some possible areas, marked 
as 1 to 4, on which subarrays could be placed. The locations are detailed in Figure 27. For 
simplicity, it is assumed that these four subarrays are all flat and rectangular in shape. 
The number of elements is 12, 6, 7 and 13, respectively (along the x-axis). The physical 
limitations in this case are significant, and the number of elements is so few that a 
straight forward combination of the four areas does not yield any reasonable performance 
at VHF. The placement of elements in each subarray needs to be optimized by the 
Minimax method in Section A to decrease the level of the grating lobes. 

By fixing the position of the edge elements in each subarray, the Minimax 
algorithm finds the optimized positions for minimum peak SEE in the array factor. Since 
subarrays 1 and 2 and subarrays 3 and 4 overlap along the x-axis, and the maximum 
possible distance is between subarrays 1 and 4, there are only two possible combinations 
of subarrays to produce a low SEE pattern for the requirement of high resolution. The 
first combination is to process all four subarrays coherently as a monostatic radar. A plot 
of the array factor and the equivalent two-way pattern is shown in the Eigures 28 and 29 
for combination 1. 

A second possible combination is to use subarrays 1 through 3 as the receive 
pattern, and then synthesize a transmit pattern for subarray 4 by the equiripple method 
discussed in the Section A. A plot of the receive pattern (3 subarrays) is given in Eigure 
30, and the two-way pattern in Eigure 31. The optimized element positions for the 2 
combinations are shown in the Tables 1 and 2. 
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Figure 26. Side view of the Aegis eruiser 
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Figure 27. 


Possible loeations of VHF subarrays on the Aegis cruiser 
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Figure 28. Combination 1: Pattern of the 4 subarrays proeessed eoherently to form a 
single array used for both transmit and reeeive. 



Figure 29. Combination 1; Equivalent 2-way pattern of the monostatic array. 
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Figure 30. Combination 2; Receive pattern using 3 subarrays 



Figure 31. Combination 2: Two-way pattern (3 receive subarrays and 1 transmit 
subarray). 
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Element position 
in A 

Subarray 1 

Subarray 2 

Subarray 3 

Subarray 4 

Xi 

0 

5.1090 

10.7439 

11.5530 

X 2 

0.0073 

6.0490 

10.7659 

11.5588 

X 3 

0.5146 

6.5390 

11.2379 

12.0647 

X 4 

2.0219 

7.0290 

11.7099 

12.5705 

^5 

2.5292 

7.5190 

12.1819 

13.9289 

^6 

3.0365 

7.6070 

12.6539 

14.5822 

Xj 

3.5438 


13.5760 

14.6629 

Xg 

4.0511 



14.9961 

X 9 

4.5071 



15.0996 

^10 

5.0657 



16.6055 

Xii 

5.5730 



17.1113 

X 12 

5.5800 



17.6171 

Xi3 




17.6230 


Table 1. Optimized loeations of elements in the 4 subarrays for eombination 1. 


Element position 
in A 

Subarray 1 

Subarray 2 

Subarray 3 

Subarray 4 

Xi 

0 

5.1090 

10.7439 

0 

X 2 

0.0073 

6.0490 

10.7659 

0.5 

X 3 

0.5146 

6.5390 

11.2379 

1 

X 4 

2.0219 

7.0290 

11.7099 

1.5 

^5 

2.0464 

7.5190 

12.1819 

2 

^6 

2.8683 

7.6070 

13.5539 

2.5 

Xj 

3.5438 


13.5760 

3 

Xg 

3.6610 



3.5 

X 9 

3.6706 



4 

^10 

4.4881 



4.5 

Xii 

5.5730 



5 

X 12 

5.5800 



5.5 

Xi3 




6 


Table 2. Optimized loeations of elements in the 4 subarrays in eombination 2. 
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Combination 1, which uses all 4 subarrays, gives better performance for both the 
peak SLL (-22.5 dB) and the beamwidth (2°) for the round trip pattern as shown on 
Figure 29. However, Combination 2 has a much lower average SLL than Combination 1, 
and the individual use of subarray 4 gives more flexibility on the design of two-way 
pattern. 

These two distributed VHF subarray designs demonstrate the possibility of 
maximizing the resolution given the spaee limitation onboard the ship. This simple 
algorithm of pattern design does not show the full potential of unequal sized or 
nonuniform spacing subarrays. Other suitable algorithms, like the Genetic Algorithm 
(GA), may provide better solutions for these types of problem [Ref 16]. 

4. Calculation of Antenna Parameters on MFAR DSA Design 

Sinee the DSA has a large spaeing between subarrays, the simple formulas for 
antenna parameters need modifieation from those for eonventional filled arrays. The 
following paragraph will show analytieal and simulated results of these important 
parameters of the MFAR DSA. The pattern parameters will be examined for the antenna 
design described in Seetion 2. 

(a) Two-way beamwidth between first nulls 

Since the result of the two-way pattern is the produet of transmit and receive 
patterns, by the prineiple of pattern multiplieation, the two-way BWFN is a funetion of 
both patterns. The BWFN of the reeeive pattern, from Equation (2), is mainly dependent 
on the subarray configuration factor, which varies as sin(A^^/2) in the uniform DSA 
design. The BWFN of transmit pattern, however, is determined by the Chebyshev 
coefficients, which are controlled by the sidelobe level, number of elements, interelement 
spacing and the frequeney. The transmit BWFN is not easy to ealeulate analytieally, 
especially for the two-way pattern. Theoretically, since the beamwidth of the receive 
pattern is sharper than that of transmit pattern, the shape of the mainlobe at broadside is 
determined mainly by the receive pattern. The BWFN of the receive pattern is 
approximately 3.2° at the H-plane broadside measured from the pattern plot on Figure 
17. 
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(b) Two-way half power beamwidth 

Since the two-way mainlobe is determined mainly by the receive pattern, the 
HPBW of the two-way pattern can be approximately determined by the receive pattern 
also, which is 




sin(^) sin(^) 


M sin(^) N sin( y) 


= 0.5, 


(14) 


It is assumed uniform weighted on both subarrays and elements. Since Ny/:^ , and 

the phase term of subarray configuration is changing much faster than the phase of 
subarray factor, the HPBW is determined mainly by the subarray configuration factor. 


A = wavelength; D = aperture wuleth 


Type of distribution, |z| < 1 

Relative 

gain 

Half-power 

beamwidth, 

deg 

Intensity of first sidelobe, 
dB below maximum intensity 

Unifomi; A(z) = 1 

1 

51A/0 

13.2 

Co.sine; /l(^) = co.s"(■77^/2): 

n = 0 

1 

5)A/D 

13.2 

n = 1 

0.810 

69 A/D 

23 

;i = 2 

0.667 

83A/D 

32 

n = 3 

0.575 

95A/D 

40 

n = 4 

0.515 

IIIA/D 

48 

Parabolic; /1(c) - 1 — (1 - A)c^: 

A = 1.0 

1 

51 A/D 

13.2 

A = 0.8 

0.994 

53A/D 

15.8 

A = 0.5 

0.970 

56A/D 

17.1 

A = 0 

0.833 

66A/D 

20.6 

Triangular; A(z) = 1 — |;| 

0.75 

13X10 

26.4 

Circular; /t(z) = V1 - 

0.865 

58.5A/D 

17.6 

Cosine-squared plus pedestal; 

0.33 + 0.66 cos^ (ttzJI) 

0.88 

63X10 

25.7 

0.08 0.92 co.s^ {'iTzJ'l), Hamming 

0.74 

76.5A/D 

42.8 


Table 3. Pattern characteristics produced by various aperture distributions [from 
Ref 10]. 


48 



For a Hamming window distribution, the HPBW is approximately 16.5X1 L from 
Table 3. For = 16, D = 5X, T = 15x5 + 2 = 77/1, by the parameters of MFAR DSA 
design in Seetion 2, the broadside HPBW is approximately 1°. This value is a little larger 
than the measured value of 0.9° from the two-way pattern on Figure 19 beeause of the 
multiplieation between the transmit and reeeive pattern. Compared to a fully populated 
uniform array (X/2), it is neeessary to have 115 linear uniform weighted elements to 
have the same resolution. It is aehieved using 80 linear elements with the DSA approaeh. 
(e) Aperture effleieney and direetivity 

Sinee the reeeive DSAs are Hamming weighted at the subarray level, the taper 
effleieney (reeeive aperture effleieney) is 0.47 using Equation (3). The effleieney of the 
transmit pattern by the same equation is 0.6. However the direetivity eannot be ealeulated 
by the equation D = 4;/r7^Acos0^ / X ^, beeause grating lobes exist. The area A is the area 
oeeupied by array elements. 

By definition, the direetivity is the ratio of the radiation intensity to the average 
radiation intensity, or [Ref. 7] 

Att 

D(0,^) = —, (15) 


l7t 7t ^ 

smOdOd^ is the beam solid angle, and E„orm is the normalized 

0 0 

eleetrie field intensity. Negleeting the element faetor, Enom is the normalized total array 
faetor. By numerieal integration of the normalized array faetor in spherieal eoordinates, 
the direetivity ealeulated from Equation (15) is about 40 dB for a uniform weighted DSA, 
whieh is almost the same value ealeulated from the same number of elements (80 x 40 ) in 
eontiguous spaeing {XU). Considering the tapering, the direetivity of the Hamming 
weighted DSA is about 36.7 dB, the direetivity of the transmit arrays (35x35 elements) 
is about 32.7 dB, and the equivalent one-way antenna gain G = ^Jg^ =^34.7dB. This 

result shows that the gain of DSA antenna depends direetly on the number of elements 
and taper effleieney. 


where 


= 11 
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Table 4 summarizes the antenna parameters for a DSA (receive subarray) design 
{N = \6, M = 5, D^=5X, Dy=5X, d^ = 0.5X, dy=0.5X). An effective area of 

(0.5/1)^ is assumed for each element. The results show that the approximate formula for 
directivity {D = ATurjAcosOJX^) is not correct in the DSA application, unless the 
occupied area is used. 



Approximate formula 

(uniform weight) 

Actual computed 

(uniform weight) 

Hamming weight 

on suharrays 

BWFN 

1.43° 

1.43° 

3.25° 

HPBW 

0.66° 

0.64° 

0.94° 

DIRECTIVITY 

45.5 dB (total area) 

40 dB (active area) 

40 dB 

36.7 dB 

MAX SLL 

-13.2 dB 

-13.2 dB (SLL) 

-3.8 dB (grating) 

-39 dB (SLL) 

-3.8 dB (grating) 


Table 4. Summary of antenna pattern parameters for MFAR DSA. 


C. SUMMARY 

This chapter has presented the basic theory and methods of subarraying. There are 
two approaches to reducing grating lobes due to the large spacing between subarrays in a 
periodic DSA. The first approach is placing the nulls of one factor ( AF^ and AFj^ ) in the 

direction of grating lobes of the other factor (AF^). The second approach is to use 

irregular spaced subarrays to reduce the peak grating lobes by redistributing the energy 
into sidelobe regions. 
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Two possible applications of DSA onboard a ship have been introduced in this 
chapter. MFAR DSA can be used at higher frequencies (X or C band) to reduce the cost 
and weight of modern radar. At lower frequencies (HF or VHF band), the application of 
DSA can increase the resolution obtained from utilizing the limited small areas on a ship. 
This is important for detecting and tracking targets with small RCS like cruise missiles. 

The calculation of antenna parameters of the DSA design shows the tradeoff in 
performance on the directivity and beamwidth. The directivity depends on the number of 
elements and tapering only; the beamwidth decreases as the subarray spacing increases. 
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IV. CONCLUSION 


A. SUMMARY 

The concept of distributed subarray antennas has been proposed for both the 
MFAR and VHP applications. By combining distributed subarrays on the available areas 
of a constrained platform, the MFAR or VHF DSA can achieve the maximum resolution 
(aside from synthetic aperture approaches) and potential reductions in cost and 
complexity. The two-way pattern design of a DSA effectively suppressed the undesired 
grating lobes by the design of separate transmit and receive patterns. From the pattern 
multiplication principle, the grating lobes in the subarray configuration pattern (AF^ ) 

have been suppressed by the design of subarray pattern ( AF ^) and transmit pattern 

The design examples of shipboard MFAR have shown that the HPBW decrease 
(increased resolution) can be achieved by spreading a fewer number of elements over a 
longer baseline and then suppressing the grating lobes with other pattern factors. The 
nulling effect of the subarray factor does not change as the beam scans away from 
broadside. The limitations and advantages of DSA have been mentioned previously 
throughout the thesis, and are summarized here. 

1. Advantages of the MFAR DSA 
a. Support of Wideband ADBF 

The collection of elements is steered in phase at the element level, and a 
beamforming network combines the element into subarrays. The subarrays are then 
steered via photonic time delay devices. This provides the necessary wideband beam 
steering and ADBF capabilities at a reasonable cost. 


53 



b. Support of Multi-functionality 

Since each subarray is physically separated, a separate function can be 
assigned each subarray more directly. The design of subarrays can be more adaptive to 
the main function from the design stage, which decreases the effects of compromises on 
performance. 

c. Low Complexity of Manufacturing and Computing with Less Control 
Elements 

Although some compromises between performance and periodic spacing 
have been made, regularity makes the implementation of subarrays more realistic and 
practical for the requirements of modern radar. However, irregular placement of elements 
is not a huge disadvantage. Randomly thinned arrays have been used for radars in the past 
[Ref. 16]. 

d.Achieving High Resolution with Less Space Limitation 
The required angular accuracy and range resolution can be achieved with 
little limitation on platform space. Little perturbations in subarray locations do not have 
much effect on the array pattern. The arrangement and dimensions of the subarrays can 
be adjusted to the shape of platform. 

e. Possibility of Multi Band Shared Apertures 

Since the spacing between the subarrays is large, it is possible to insert 
other low frequency elements to share the same aperture space. For example, if the 
subarrays structure is designed for X-band frequencies, then an L-band array element can 
be inserted between subarrays with half wavelength spacing. Figure 32 shows a possible 
scheme for a X and L band array. 
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X band subarray 


L band 
element 


L band 
element 


X band subarray 



PHASE 

SHIFTERS 

SUBARRAYS 
IN/OUT 


LOW NOISE 
AMPLIFIER 


TRUE 

TIME 

DELAY 


ARRAY IN/OUT 


L/X BAND 


Figure 32. Possible beam eontrol seheme of a dual-band DSA [after Ref 17] 


2, Limitations of the MFAR DSA Design 

a. Less Efficiency on Surveillance Due to the Narrow Beamwidth 
of Transmit Pattern 

If the beamwidth of the transmit pattern is tied to the synthesis of the two- 
way pattern, the resulting beam may not provide a large enough field of view for 
searehing. This disadvantage, however, can be corrected by dividing the transmit pattern 
into searching and tracking modes. In the search mode, the center transmit/receive 
module operates alone like any other element-based phased array. In the tracking mode, 
the center transmit module is a high contrast pencil beam, and the resultant two-way 
pattern is the multiplication of transmit and receive pattern synthesized for the DSA. A 
dual-band shared aperture is another possible approach for separating the track and 
surveillance into different combinations; for examples, long-range surveillance (L band), 
mid-range surveillance (X band) and tracking (X band DSA). 
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b. Lower Beam Efficiency Due to the Grating Lobes 

From the calculations in Chapter III, the antenna directivity depends only on 
the area occupied by elements. There is a reduction in beam efficiency resulting from the 
grating lobes, and any sidelobe taper will lower the efficiency further (typically about 3 
dB). However, this can be improved at the expense of more subarrays or closer spacing 
between them 

c. Fewer Degrees of Freedom on Pattern Synthesis 

Since the number of individual control elements is decreased by the subarray 
design, the degrees of freedom (DOF) in terms of ADBF capability is decreased 
substantially. The limitations resulting from the synthesis of transmit and receive patterns 
also restrict the implementation of DSA in the environment of high clutter because of 
lower sensitivity and less capability for nulling interferences. 

d. More Complexity of Task Scheduling 

Having a multifunction radar complicates the task scheduling of a single 
antenna. The two-way pattern approach introduces more constraints on the scheduling 
tasks. The time budget i.e., the allocation of radar time to different tasks is dependent on 
the radar antenna parameters. Since the main tasks of MFAR are search, track and other 
auxiliary functions, the strategy of how to execute of all tasks in the best possible way is 
very complicated. 


B, POSSIBLE FURTHER RESEARCH TOPICS 

There are two primary directions of research on suppressing the grating lobes in 
the pattern of widely spaced subarrays, as discussed below. 

1. Filtering Approach 

There are many filtering techniques already developed in the area of digital signal 
processing. If the main objective is to suppress the grating lobes from self-interference, 
the adaptive notch filter (ANF) or space-time domain adaptive processing (STAP) might 
be the possible research area for the implementation of DSA. 
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Adaptive notch filters are widely used in many signal-processing applications to 
extract and trace the narrow-band noise. The basic principle is to produce the adaptable 
frequency response in both the time and space domain: zero at the specified spatial angle, 
one otherwise. Figure 33 shows the example of frequency response on nulling the angles 
±15° and ±37°, everywhere else unchanged. The depth and width of nulls should be 
adjusted according the grating conditions. But this must be implemented in both the space 
and time domain, which introduces the applications of the STAP algorithm [Ref. 18]. For 
a DSA the notches would be placed at grating lobe locations. 

STAP processes signals in the spatial and time domains. The time domain 
includes both slow-time (pulse repetition interval) and fast-time (range cell). Therefore, 
the utilization of phase and amplitude weighting as a means to achieve a desired steering 
or nulling direction must consider the spatial, slow-time and fast-time factors. 

2, Digital Arrays 

Unlike an analog beamforming network, digital beamforming arrays digitize 
received signals at the element level, then process these signals in a digital processor to 
form the desired beam and frequency response (in this case, nulling in the grating 
direction). By doing so, the total information available at the aperture has been properly 
preserved and can be manipulated indefinitely without introducing further error (other 
than computer roundoff error). 

Adaptive digital beamforming at the element level can reject interference and at 
the same time steer a main lobe in the direction of a desired signal. Therefore, any 
number of beams can be formed or the signal can be rejected according to a selected 
algorithm. 
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Figure 33. Adaptive nulling at the angle ±15°, ±37° by ANF. 
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APPENDIX A: GLOSSARY OF TERMINOLOGY 


MFAR 

ADBF 

LPI 

MFSAR 

dB 

mrad 

MHz 

GHz 

MESAR 

TWS 

AF 

AF^ 

AF^ 

AFr. 

AF^ 

EF 

BWFN 

HPBW 

A 

SNR 

RCS 


Multi-function Array Radar 
Adaptive Digital Beamforming 
Eow Probability of Intereept 
Multifunetion Subarray Radar 
Deeibels 

milli radian (10“^ radian) 

Megahertz (10*’ cyeles/seeond) 

Gigahertz (10‘’ eyeles/seeond) 

Multifunetion Eleotronieally Seanned Adaptive Radar 

Traek-While-Sean 

Array Factor 

Subarray Factor 

Subarray Configuration Factor 
Transmit Array Factor 
Receive Array Faetor 
Element Eaetor 

Beamwidth Between Eirst Nulls 
Half Power Beamwidth 
Eambda — wavelength 
Signal-to-Noise Ratio 
Radar Cross Seetion 
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SLL 

Sidelobe Level 

DSA 

Distributed Subarray Antenna 

FIR 

Finite Impulse Response 

WLS 

Weighted Least Squares algorithm 

Tx 

Transmitter 

Rx 

Reeeiver 

HF 

High Frequeney (3 ~ 30 MHz) 

VHF 

Very High Frequeney (30-300 MHz) 

UHF 

Ultra High Frequeney (300 - 3000 MHz) 

DOF 

Degree Of Freedom 

ANF 

Adaptive Noteh Filter 

STAP 

Spaee Time Adaptive Proeessing 
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APPENDIX B: MATLAB CODE LISTING 


The major MATLAB Code listing for the pattern plot as the following. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% DSA.m 

% Pattern calculation and plot for Figure 8,13,14,18-25 
% Quasi-monostatic subarrays AF 

% XMTR AF has 35x35 Remez-weighted elements spacing 0.5 wavelength 
% RCVR is 16x8 subarrays with 5x5 elements in each subarray 
% Each subarray spacing 5 wavelengths in x&y axis with Flamming weights 
% Display the aperture efficiency also 
% resolution of plot is quarter degree 

clear all; 

llx=0:5:75; % 16 Subarray spaced (Lambda) in x-axis 
lly=0:5:35; % 8 Subarray spaced in y-axis 

dx=0.5; Mx=5; % element spacing and number in subarray along x-axis 
dy=0.5; My=5; % element spacing and number in subarray along y-axis 

% Scan angle input 

% Figure 21-25 need to change angle accordingly 

thets=0; 

phis=0; 

us=sin(thets*pi/180)*cos(phis*pi/180); 

vs=sin(thets*pi/180)*sin(phis*pi/180); 

thet=hnspace(-90,90,720); % x,y direction cosine 
theta=thet.'; 

phi=hnspace(0,180,720); 
u=sin(theta*pi/l 80)*cos(phi*pi/l 80); 
v=sin(theta*pi/180)*sin(phi*pi/180); 

U=u-us; 

V=v-vs; 

% Subarray configuration AF 
% Weighting Function by Flamming window 
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wx=zeros(l,length(llx)); % weights in x-axis 
for n=l :length(llx) 

xn=(n-1 )/(length(llx)-1 )*2-1; 
wx(n)=0.92*abs(cos(xn*pi/2))''2+0.08; 
end 

wy=zeros(l,length(lly)); % weights iny-axis 
for n=l :length(lly) 

yn=(n-1 )/(length(lly)-1 )*2-1; 
wy(n)=0.92 *abs(cos(yn*pi/2))''2+0.0 8; 
end 

CAP x=zeros(720); 
for lx=l :length(llx) 
psilx=2 *pi *llx(lx) *U; 

CAP x=C AP x+wx(lx) *exp(j*psilx); 
end 

CAPy=zeros(720); 
for ly=l :length(lly) 
psily=2*pi*lly(ly)*V; 
CAPy=CAPy+wy(ly)*exp(j *psily); 
end 

CAP=CAPx.*CAPy/sum(wx)/sum(wy); 

% Single subarray AP 
SAPx=zeros(720); 
for sx=l:Mx 
psix=2*pi*dx*U; 
SAPx=SAPx+exp(j*(sx-l)*psix); 
end 

S APy=zeros(720); 
for sy=l:My 
psiy=2*pi*dy*V; 

S APy=S APy+exp(j * (sy-1) *psiy); 
end 

SAP^SAPx. * S APy/Mx/My; 
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% XTMR array design to pass mainlobe, others with 45 db attenuation 
% First ealeulation the equal ripple eoeffieients 

% use for obtaining 9 deg mainlobes, every other angle response -40 db 

px=35; % Element numbers in x-axis 

fs=90; % Array faetor half spaee 

rp= 1; % Passband ripple in db 

rs=-50; % Stopband desired attenuation in db 

fe=[0 4.5]; % Cutoff spatial freq 

a=[l 0]; % Desired amplitude 

dev=[(10^(rp/20)-l)/(10^(rp/20)-l-l) 10^(rs/20)]; % dB to numerieal 
[nx, f0x,a0x, W x]=remezord(fe,a,dev, fs); 
bx=remez(px-l,fOx,aOx,Wx); % Using fixed number elements 
bx=abs(bx(l:length(bx))); % Elements weighting 
Bx=bx/max(bx); % Normalized weights 

py=35; % element # in y-axis 

fs=90; % Array faetor half spaee 

rp=l; % Passband ripple in db 

rs=-50; % Stopband attenuation in db 

fC=[0 4.5]; % Cutoff spatial freq(angle) 

a=[l 0]; % Desired amplitude 

dev=[( 10^(rp/20)-1)/(10^(rp/20)-^ 1) 10^(rs/20)]; 

[ny,fOy,aOy,Wy]=remezord(fC,a,dev,fs); 

by=remez(py-l,fOy,aOy,Wy); % Using fixed number elements 

by=abs(by(l:length(by))); % elements weighting 

By=by/max(by); % Normalized weights 

% Array faetor in square XTMR 
XAFx=zeros(720); 
for Xx=l :length(Bx) 
psiXx=2*pi*dx*U; 

XAFx=XAFx-KBx(Xx)*expG*(Xx-l)*psiXx); 

end 

XAFy=zeros(720); 
for Xy=l :length(By) 
psiXy=2*pi*dy*V; 
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XAFy=XAFy+By(Xy)*expG*(Xy-l)*psiXy); 

end 

% calculate RCVR aperture efficiency 
Rsl=0; Rs2=0; 
for il=l:length(llx) 
for i2=l :length(lly) 

Rs 1 =Rs 1 +abs(wx(i 1)*wy(i2)); Rs2=Rs2+abs(wx(il)*wy(i2))^2; 
end 
end 

eta=Rs 1 ^2/length(llx)/length(lly)/Rs2; 

% calculate XTMR aperture efficiency 
Xsl=0; Xs2=0; 
for il=l:length(Bx) 
for i2=l :length(By) 

Xsl=Xsl+abs(Bx(il)*By(i2));Xs2=Xs2+abs(Bx(il)*By(i2))"'2; 

end 

end 

ETA=Xs 1 ^2/length(Bx)/length(By)/Xs2; 
disp(['RCVR aperture efficiency: ',num2str(eta)]) 
disp(['XTMR aperture efficiency: ',num2str(ETA)]) 

% XTMR AF 

XAF=(XAFx. *XAFy)/sum(Bx)/sum(By); 
dbXAF=20*logl0(abs(XAF)); 

% RCVR AF 

RAF=SAF.*CAF; 

dbRAF=20*logl0(abs(RAF)); 

% Two-way pattern 
BAF=XAF.*RAF; 

dbBAF=20*logl0(abs(BAF)); % 2-way pattern in dB 
dbSAF=20*logl0(abs(SAF)); % Subarray pattern in dB 
dbCAF=20*logl0(abs(CAF)); % Subarray configuration pattern in dB 

figure(l) % Figure 8 

subplot(3,l,l), plot(theta(l:720), dbRAF((l:720),phis-M)),axis([-90 90 -60 0]), 
grid on, title(' Receive pattern of 16 subarrays spacing 5\lambda') 
subplot(3,l,2), plot(theta(l:720), dbSAF((l:720),phis-M)),axis([-90 90 -60 0]), 
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grid on, title(' Pattern of single subarray which is composed by 5 elements') 
ylabel('Relative pattern at \phi=0\circ in dB') 

subplot(3,l,3), plot(theta(l:720), dbCAF((l:720),phis+l)),axis([-90 90 -60 0]), 
grid on, title(' Pattern of 16 subarray configuration') 
xlabel('\theta in degree') 

figure(2) 

subplot(3,l,l), plot(theta(l:720), dbBAF((l:720),phis-H)),axis([-90 90 -100 0]), 
grid on, title(' Two Way pattern of 16 subarrays spacing 5\lambda') 
subplot(3,l,2), plot(theta(l:720), dbRAF((l:720),phis-H)),axis([-90 90 -100 0]), 
grid on, title(' Receive pattern of 16 subarrays spacing 5\lambda') 
ylabel('Relative pattern at \phi=0\circ in dB') 

subplot(3,l,3), plot(theta(l:720), dbXAF((l:720),phis-H)),axis([-90 90 -100 0]), 
grid on, title(' Transmit pattern of 35 elements') 
xlabel('\theta in degree') 

figure(3) % Figure 18 

% Receive pattern of 16 subarrays spaced 5 wavelength 
plot(theta(l:720), dbRAF(l:720,phis-H)),axis([-90 90 -60 0]), 
grid on, 

ylabel('Receive pattern at \phi = 0 deg in dB') 
xlabel('\theta in degrees') 

figure(4) % Figure 13 
% Transmit pattern of 35 elements 

plot(theta(l:720), dbXAF(l:720,phis-H)),axis([-90 90 -80 0]), 
grid on, 

xlabel('\theta in degrees') 

ylabel('Transmit pattern at \phi = 0 deg in dB') 

figure(5) % Figure 19 

% plot Two-way pattern of 16x8 subarrays spaced 5 wavelength at Fl-plane 
plot(theta(l:720), dbBAF(l:720,l)),axis([-90 90 -80 0]), 
grid on,ylabel('Two-way pattern at \phi = 0 deg in dB'), 
xlabel('\theta in degrees') 

figure(6) % Figure 20 
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% plot Two-way pattern of 16x8 subarrays spaced 5 wavelength at E-plane 
plot(theta(l:720), dbBAF( 1:720,36l)),axis([-90 90 -80 0]), 
grid on,ylabel('Two-way pattern at \phi = 90 deg in dB'), 
xlabel('\theta in degrees') 

% End of DSA.m 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% IrregSubarray.m 

% AF of VHF DSA with freq 100 MHz 

% Looking for the locations of 38 optimized element in 4 subarrays 
% Find element optimized location and first null by myfunlO.m 
% All position in wavelengths 
% Used for Figure 28, 29 pattern plot 

clear all; 

u0=0.05; % Initial guess of beamwidth 
ul=l; 

% element locations of first subarray 

s0=0; % Lower boundary of elements in 1st subarray 

si 1=5.58; % Upper boundary of elements 

for n=l:10 

sl(n)=0.5073*n; 

end 

% Initial guess and boundary of 12 element locations 
sl0=cat(2,s0,sl,sl 1); 
sll=cat(2,s0,sl-0.5,sl 1); 
slu=cat(2,s0,sU-0.5,sl 1); 
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% element loeations of 2nd subarray 
s21=5.109; % First element position 
s26=7.607; % Last element position 
for n=l:4 

s2(n)=s21+n*0.49; 

end 

% Initial guess and boundary of 6 element loeations 
s20=eat(2,s21 ,s2,s26); 
s21=eat(2,s21,s2-0.45,s26); 
s2u=cat(2,s21,s2+0.45,s26); 

% element loeations of 3rd subarray 
s31=10.7439; % First element position 
s37=13.576; % Last element position 
for n=l:5 

s3(n)=s31+n*0.472; 

end 

% Initial guess and boundary of 7 element loeations 
s30=eat(2,s31,s3,s37); 
s31=eat(2,s31,s3-0.45,s37); 
s3u=cat(2,s31,s3+0.45,s37); 

% element loeations of 4th subarray 
s41=11.553; % First element position 
s413=17.623; % Last element position 
for n=l:ll 

s4(n)=s41+n*0.50583; 

end 

% Initial guess and boundary of 13 element loeations 
s40=eat(2,s41,s4,s413); 
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s41=cat(2, s41, s4-0.5, s413); 
s4u=cat(2,s41,s4+0.5,s413); 

% Combination of all element positions 

Ilb=eat(2,u0,sll,s21,s31,s41); % lower bound of element and the first null position 
Iub=eat(2,ul,slu,s2u,s3u,s4u); % upper bound of element and the first null position 
10=eat(2,u0,sl0,s20,s30,s40); % Intial guess of element and the first null 
options=optimset(']VlinAbsMax', 1); 

[l,fval, maxfval, exitflag, output]=fminimax(@myfunlO, 10 ,[],[];[];[] 41 b,lub,[],options) 

% Result pattern ealeulation 
% Sean angle input in degrees 
thets=0; phis=0; 

us=sin(thets*pi/180)*eos(phis*pi/180); 

vs=sin(thets*pi/180)*sin(phis*pi/180); 

% x,y direetion eosine 

thet=hnspaee(-90,89.75,720); 

theta=thet.'; 

phi=hnspaee(0,179.75,720); 
u=sin(theta*pi/l 80)*eos(phi*pi/l 80); 
v=sin(theta*pi/180)*sin(phi*pi/180); 

U=u-us; V=v-vs; 

% AF of optimized loeation 
LAFx=zeros(720); 
for lx=l:length(ll) 
psilx=2 *pi *ll(lx) *U; 

LAFx=LAFx+exp(j *psilx); 
end 

dbLAF=20*logl0(abs(LAFx)/length(ll)); 
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db2LF=40*logl0(abs(LAFx)/length(ll)); % Equivalent 2way pattern 


% Pattern plot of optimized position used on Fig. 28 
figure(l) 

plot(theta, dbLAF(:,4*phis+l)),axis([-90 90 -60 0]), 
grid on, 

xlabel('\theta in degrees') 
ylabel('Relative pattern in dB') 

% Equvalent 2-way pattern used on Fig. 29 
figure(2) 

plot(theta, db2LF(:,4*phis-l-l)),axis([-90 90 -60 0]), 
grid on, 

xlabel('\theta in degrees') 
ylabel('Relative pattern in dB') 

% end of IrregSubarray.m 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% myfunlO.m 

% Function for optimization location of subarrays 
% Function called by the IrregSubarray.m 

% Minimaxing the peak sidelobe level in the visible region[-pi/2,pi/2] 

function F=myfunl0(l) 

% Array factor of same subarray with irregular spacing 
% Computation for U axis only 
u0=linspace(l(l),0.5*pi,5000); 


alpha=2*l; 
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F1 =zeros( 1 ,length(uO)); 
for 11=2:39 

F1 =F 1 +exp(j *alpha(n) *u0); 
end 

F=abs(Fl); 

% end of myfunlO.m 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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